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Abstract 

Thermal models have proven to be an useful and simple tool used to make theoretical predictions 
and data analysis in relativistic and ultra-relativistic heavy ion collisions. A new version of these 
models is presented here, incorporating a non equilibrium feature to the description of the inter- 
mediate fireball state formed at the chemical freeze-out. Two different effective temperatures are 
attributed to the expanding fireball, regarding its baryonic and mesonic sectors. The proposal is 
not merely to include an additional degree of freedom to reach a better adjustment to the data, but 
to open a room in the model conception for considerations on the non-equilibrium scenario of the 
system evolution. A set of well consolidated data for particles production is used to validated the 
reformulated version of thermal models presented here. A rather good performance of the extended 
version was verified, both for the quality of particle ratio data fittings as well as for describing the 
asymptotic energy behavior of temperatures and baryochemical potential of the colliding nuclear 
system. 

PACS numbers: 25.75.-q,13.85.-t,13.75.-n 
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I. INTRODUCTION 



A significant amount of experimental data for particle production in liigfi energy fieavy 
ion collisions (10 - 200 GeV / A at center of mass) has been accumulated during the last 
two decades. Many different theoretical attempts have tried to described these data using 
thermal models in the approximation of global thermal equilibrium, considering just one 
freeze-out temperature P-[16]. Because of its simplicity, these models are widely used to 
treat and explain the data. However, thermal models often are not able to describe ade- 
quately all multiplicities of hadrons [4J . For instance, the abundance of strange particles are 
overestimate and the pion yields are underestimated. In addition to that, particle number 
density and total energy density calculated in chemical freezing are overestimated as well 
[To] . To overcome these limitations some features have been incorporated into the model, 
like the inclusion of a phenomenological parameter 7^ to adjust the strangeness production 
in the reaction process [5]-[7], or even the inclusion of effects of excluded volume p]-|ll]. 
However, the main criticism is that the oversimplified conception of the model makes its con- 
nection with more elaborated models [I7j-|22] a hard task. This happens for example when 
trying to compare results from thermal models with those from hydrodynamical models, 
non equilibrium kinetic equation approach or even with others phenomenological statistical 
model including the formation of separated thermal subsystem [8]. Non-equilibrium pro- 
posals for a phenomenological treatment of the hadronic system has been recalled with the 
advent of lattice calculation results, showing that the hadronization critical temperature 
differs significantly from the chemical freeze-out temperature estimated, using thermal con- 
ventional models, by tens of MeV [21]- [26]. The inclusion of non-equilibrium features in 
handling with hadron production in high energy collision were revived trying to justify the 
difference between critical and freeze out temperatures in these calculations. However, the 
used arguments were not strong enough to be compatible with RHIC results few years latter 
[27]- [SU]. The solution is still far from being conclusive, but some effort have been made 
more recently in this direction [31J. The main idea goes back to the Hagedorn's states, 
introduced in the 1960's, to quickly drive hadrons to chemical equilibrium [27J. The intro- 
duced gas of hadronic resonant states should have a mass spectrum parametrized by the 
Hagedorn's exponential form using a critical temperature, the Hagedorn temperature. The 
discussion of the possibility of have different values of Hagedorn temperature for baryons 
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and mesons in the hadronization process was addressed in Ref.[28]. However, the inclusion 
of the highest mass branch of observed hadronic resonant states in the Hagedorn exponential 
parametrization of the spectrum fails, and the critical temperature parameter value can not 
be well defined |32]-[33]. 

In the present work we assume that the chemical freeze-out stage is not characterized by 
a global thermal equilibrium with a single temperature of the expanding fireball. Instead 
of this, it was introduced a picture of the fireball freeze-out stage out of equilibrium hav- 
ing baryonic and mesonic sectors characterised by two distinct temperatures, which could 
be considered as being an effective thermodynamic quantities associated to these sectors. 
Indeed, we have no well founded reason to support the conventional thermal model state- 
ment of equal thermal evolution for mesons and baryons during the whole evolution of the 
system. The chemical freeze-out can occur in non-equilibrium condition before the kinetic 
freeze-out regime. Only after the later stage it is believed that the thermodynamic equilib- 
rium is reached. Consequently, only simplicity arguments supports the isothermal fireball 
assumption for the chemical freeze-out. 

Since the pioneering works by Fermi [35] and Heisenberg [36] on thermal models, the idea 
of treating baryons and mesons as distinct thermal sectors was suggested. In these works 
baryons arise merged in a meson cloud and the description for particle multiplicity are 
essentially statistical in nature. When the concept of temperature is introduced in Fermi's 
work, it is defined by relating the number of particles with the energy density of the system 
by mean of the Stefans-Boltzmann law. The Hagedorn resonance mass spectra mentioned 
previously is nothing more than another way to associate an effective temperature to the 
thermal state of the system. In this case the assumption is that baryons and mesons are 
produced in thermal contact with a correspondent Hagedorn thermal reservoir, defined by 
the observable mass spectra of baryons and mesons. Different effective temperature values 
for mesons and baryons are extracted from this association [32]-[33|. In addition, due to 
the experimental difficulties in establish the observable particle resonance spectra, some 
ambiguities values of these effective temperatures can be found [33] . 

In our case, we are assuming that baryons and mesons are directly in contact with the 
thermal vacuum from where particles are produced. This is an image compatible with 
the use of a grand canonical ensemble behind the Fermi distribution (for baryons) and 
Bose-Einstein (for meson) employed to determine the density of particles at the freeze-out 
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in thermal models. At the beginning of the chemical freeze out mesons and baryons are 
produced from an initial thermal vacuum state, whereas at the final stage, and due to the 
depletion of energy density of the expanding system, the vacuum should be considered at a 
different thermal state. The meson freeze-out is associated to this final stage. Wc remark 
that these two effective temperatures introduced here in the thermal model still preserves the 
relatively simple model skill to perform calculation, bringing this kind of approach closer 
to a more detailed description of the fireball expansion. Having in mind this qualitative 
arguments, a necessary condition to validate this alternative model is the improving of 
the results obtained with the conventional calculation. Thus, as a preliminary task, we 
applied the two-tcmpcraturcs model to different collisional processes with well established 
data set for particle ratios, and comparing results with those obtained with the conventional 
model. In addition, we discuss the behavior of freeze out temperatures, entropy and baryonic 
chemical potential for increasing collision energy of nucleus-nucleus reaction. 

This work is organized presenting firstly in Section II, the main changes introduced in 
the conventional thermal model (IT-Model) due to the two-temperatures description {2T- 
Model). The results obtained for particle ratios within the modified model and they com- 
parison with the results of a conventional thermal calculation is shown in Section III. The 
analysis of temperature, entropy and chemical potential behavior with the increasing colli- 
sion energy is presented in Section IV. And finally, our conclusions and final discussions are 
given in Section V. 

II. THERMAL MODEL WITH TWO FREEZE-OUT TEMPERATURES 

In thermal models the particle density of the A;th particle is obtained using the distribution 
functions of baryons and mesons given by the Fermi- Dirac and Bose- Einstein distributions 
functions respectively, in the following way 



where = yp^ -|- m|, the minus (plus) sign stands for the Fermi-Dirac (Bose-Einstein) 

distribution, 7^ is the spin degeneracy and is the chemical potential, using Tj {i — 1,2; 
baryon,meson) for temperature of different sectors. The chemical potential of the kih. particle 
/Xfe is calculated using the global chemical potentials {//b, /xj3, ^s} ^ IJ'k — B^/ib + I^kl^'is + 




(1) 
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SkfJ'S, where {Bk, I3k, Sk} are the baryonic, isospin and strangeness numbers of the particle 
respectively. The global chemical potentials are constrained to the conservation of global 
baryonic, isospin and strangeness charges as 

VY^BkUk = Zi + Z2 + Ni + N2 (2) 

k 

vY^n^k = \{Zi + Z2-N,- N2) (3) 
k ^ 

Y.Sknk = ^, (4) 

k 

where Zi and A^j ii = 1, 2; target, projectile) are the number of protons and neutrons of the 
colliding nuclei and V is the volume of the fireball. This procedure can be understood in 
connection with the Generalized Gibbs Ensemble prescription [31]. 

The complete lightest baryonic octet, decuplet, and mesonic nonet are included, together 
with around two hundred of other baryonic and mesonic resonances with masses up to 2 
GeV. It is also important to point out that the particle population calculated by Eqs. ((iQ 
(to be compared to the data) should be modified due to the decay chain of resonances, a 
process known as the feed-down of particle populations. We reformulated the resonance 
decay feeding processes of the conventional thermal model in order to make it consistent 
with the two-temperatures approach. For that, the final density of produced particles are 
given by. 



n 



k 



nk + Yl ^J^k rij, (5) 



where is the thermal particle density calculated using Eq. ([T]), and Tj^i^ is the probability 
of particle j to decay into particle k, taken from experimental particles decay data. Popu- 
lations of thermal particle Uk are determined consistently with Eqs. (|2|-(|4] using in Eq. ([T]) 
the specific effective temperature for baryons and mesons. 



III. PARTICLES RATIOS RESULTS IN 2r-MODEL 

The results obtained from the conventional model with one temperature, referred in 
this work as IT-Model is now compared with those from the two-temperatures model, the 
introduced 2T-Model. We have performed the conventional calculation by using the 2T- 
Model code with the constraint Tf, = T^. The two-temperatures code was written from 
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a previous version of thermal model developed by two of the authors [T5]-[T6]. A Monte 
Carlo sampling of the model parameter values covering a wide range of the domain was 
taken in order to speed up the computational calculation. Results of the parameter values 
(temperatures and baryochemical potentia) for the best fitting of the particle population 
ratios produced in the five systems under analysis are compared in Tables [l] (IT- Model) and 
[n] (2-Tmodel). The CI data set was extracted from Ref. [12], C2 data was taken from Ref. 
[3] and C3, C4 and C5 were taken from |37]. To compare the quality of the fittings we also 
show the correspondent Xd.o.f values of the data adjustment. 

TABLE I. Results for the thermodynamic parameters for different collisions using conventional 
IT-Model. The error bars in the temperature and baryonic chemical potential are established as 
the isolated variation in these quantities which leads to a 10% of change in the optimal Xdo f value. 



Sets 


T (MeV) 


f^b (MeV) 


Xd.o.f 


^SNN [OeV) 


Ions 


CI 




779 qC^-^S) 
' '^-^(0.25) 


0.160 


2.32 


Ni + Ni 


C2 


147 4^-^? 

-^(2.6) 


907 9(14.9) 


1.667 


17.3 


Pb + Pb 


C3 


irn .(9.47) 
^^^•^(8.70) 




7.435 


62.4 


Au + Au 


C4 


156.50[?:f,] 




5.876 


130 


Au + Au 


C5 


i69.5o[;:S] 


^"^•^^(6.80) 


3.159 


200 


Au + Au 



TABLE II. Thermodynamic parameters for different collisions obtained by using the 2T-Model. 
The error bars are defined as previously in Table |l| 



Sets 


Tm {MeV) 


n {MeV) 


^lh {MeV) 


X^d.o.f 


^SNN (Cel/) 


Ions 


CI 


27 9^2-^^) 

^'•^(0.09) 


^^■"(0.48) 


779 n^^-'') 

''^•^(10.35) 


0.9666 


2.32 


Ni + Ni 


C2 


^•^"^••^"(0.44) 


140.05[?:-) 


9C1 on(6.60) 
ZOi.OU(.^2.80) 


0.9997 


17.3 


Pb + Pb 


C3 




145.50[?:Sj 


^^•^(6.35) 


0.9996 


62.4 


Au + Au 


C4 


129 74^^-^^) 

^^^•'^(0.19) 


146 32(°-^^) 


24 91^^-^^^ 

^^•^■^(0.16) 


0.9999 


130 


Au + Au 


C5 


l^S-Olloigj 


147 36^°-^^^ 

^*'-'^"(0.48) 


90 91 (0.89) 
^■^•^^(9.41) 


1.0 


200 


Au + Au 



Comparing the results we see that in all cases both effective temperatures in Table |TT] (Tf, 
and Tm) are lower than the single freeze out temperature in Table [l] (T). Another direct 
observation is that the fits in Table |TT] present a better quality when compared with those 
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in Table [l| which is reflected in the Xd.o.f values. 

In Figure [l] we show plots of the parameter values for data adjustment with chi'^ o / 
a narrow range ( 0.9 < Xd.o.f < ^■^) ^^e systems Pb+Pb at 17.3 GeV and Au+Au at 
200 GeV. Each point corresponds to values of {Tb,Tm, fib} in the parametric space. In left 
panels the difference 6T = — Tm are plotted against the baryonic chemical potential. A 
set of five hundred points are generated covering the Xdo f window. A Monte Carlo sort of 
the parameters values was performed exploring a wide domain of the parameters in order to 
speed up the computational calculation. In central panels we depicted the same difference 
against the baryonic effective temperature, and in right panel we have the point spread 
out in the (Tf,, Tm) plane. Similar plots were obtained for the other analyzed system (not 
shown) with the equivalent density of points in this narrow Xdo f window. The values of 6T 
in these plots does not reach null values {{Tb = Tm = T), indicating that IT-Model is not 
able reproduce data when the requested quality of adjustment is established. Except for 
the lightest system with lowest collision energy it was possible to have the data adjustment 
with Xd.o.f established range with both models. 

In Figures |2] and |3] we show the obtained best fitting of particles ratios results correspond- 
ing to the system in previous figure compared with the data. Although the fitting seems 
similar in the plots the rather good improvement in the quality by using 2r-Model can be 
seen comparing the Xd.o.f values in previous tables. 
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FIG. 1. Plots for thermal parameters of Pb + Pb collisions at ^snn = 17.3 GeV in panel (a). In 
panel (b) we have the same plots for Au + Au at ^/snn = 200.0 GcV. In both calculation it was 
requested that the parameters values leads to a particle data rate fitting with 0.9 < Xdof < 
In the left column we have ST vs plots, in the central column we have ST vs T^ plots and in the 
right column we show the points in the {Ti,,TTn) plane. 



9 




TI+/7I- K+zX" A/A K /7I- A/K+ A/K" 




7r+/7i- K+/K" A/A K /%- A/K+ A/K" 



FIG. 2. Plots for particle ratios of Pb + Pb at ^/s^=17.3 GeV fitted using the IT-Model, plot 
(a), and the 2T-Model, plot (b), both with open circles. The full circles with error bars are the 
data taken from Refs.[12J and [3]. 
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FIG. 3. Plots for particle ratios of Au + Au at 200 GeV fitted using the IT-Model in part (a), and 
the 2r- Model, in (b), both with open circles. Full circles with error bars are the data taken from 
Ref.[37j. 
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IV. THE COLLISION ENERGY DEPENDENCE OF THE FREEZE-OUT TEM- 
PERATURES, BARYOCHEMICAL POTENTIAL AND ENTROPY 



From the results presented in Tables |T] and [IT] we can exhibit a sketch of the baryonic and 
mesonic freeze-out temperatures and baryochemical potential behavior as a function of the 
collision energy ^/S~NN. With the five points, corresponding to the chosen systems, we fitted 
the results of the thermal parameters as function of the collision energy, a/S^vtv- For the 
temperatures we have used the fitting expression similar to that one considered in Ref . [13] , 

164c 

T = (6) 
where v = ^/Snn- For baryonic potential the parametrization used was [13] 

= (7) 

1 + ev 



The parameters values for the best fitting in Eqs. ^ and ([T]) are showed in Tables III 
and HV] for both IT and 2T models. 

In Figs. |4] and [5] it is shown a plot of temperature and baryonic chemical potential 
dependence with the collision energy. As observed in many others calculations, a sudden 
increase is observed in the temperatures (Fig. |4]) in the low energy regime and a constant 
value reached asymptotically. Our result shows a quite defined plateau for the 2T-Model 
reached quickly with the energy increasing, that is not the case for the IT-Model result in 
which it is not formed for the covered collision energy values. 

In Fig. [5] it is shown the freeze out baryochemical potential as a function of the collision 
energy for the case of IT and 2T models. As the meson sector does not carries baryonic 
charge, the results for the baryochemical potential cannot differ significantly when using (IT) 
or (2T) description. A small difference is observed due to the effect of the temperatures on 
the constraints coupling (Eqs. ([l]|4])) and on the feeding process in Eq. 

In 2T-Model the entropy for mesonic or baryonic sectors are calculated separately as. 



Si = Y (^ei + Pi-^ /ifc^fe j 



(8) 



where i = m(meson), 6(baryon) specifies the sector, with Si and Pi being the corresponding 
energy density and pressure. The k subscript refers to the particles entering in the sector 
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TABLE III. Parameters of the temperature fitting in Eq. Q. First line corresponds to the fit 
within the IT-Model, the other two lines are the fittings within the 2r-Model. 





a 


b 


c 


T 


19.026 


1.739 


1.0 


n 


24.108 


2.145 


0.897 


T 


1115.926 


6.783 


0.800 



TABLE IV. Parameters of the chemical potential fitting in Eq. Q. First line corresponds to the 
fit within the IT-Model, the second line is the fittings within the 2r-Model. 





d 


e 




1094.050 


0.178 




1177.078 


0.225 



composition. It is well known from conventional models that the dimensionless quantity 
S/T'^ increases (decreases) for mesons (baryons) as y/S^N increases This behavior is 
claimed as a change from a baryon dominance to a meson dominance regime of the collision 
freeze-out, marked by the energy value at the point where Sm/T^ and Sh/T^ curves crosses. 



From Table [VT] this change of regime is also observed in the 2T- Model, where can be identified 
a meson dominance regime around \/Snn = 7 GeV. 



TABLE V. Mesonic and baryonic entropy densities for the IT-Model. The third and fourth columns 
are the mesonic and baryonic entropy densities for different collisional sets respectively. In the fifth 
column we can see the mesonic-baryonic entropy ratio. In the last two columns we have the mesonic 
(baryonic) entropy density over cubic mesonic (baryonic) temperatures. 



Set 


y/sNN{OeV) 


Sm (fm-3) 


Sb (fm-3) 


Sm/ Sf) 


Sm/T"^ 


Sb/T^ 


CI 


2.32 


0.669 X 10-^ 


0.492 X 10^2 


0.136 


0.198 


0.146 


C2 


17.3 


1.940 


0.839 


2.310 


4.652 


2.013 


C3 


62.4 


1.994 


0.445 


4.486 


4.495 


1.002 


C4 


130 


2.368 


0.581 


4.077 


4.747 


1.164 


C5 


200 


3.230 


1.008 


3.205 


5.096 


1.590 
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FIG. 4. Plots of temperature against the collision energy. The plot with full-line and full-circles 
represents the parametrization and fitted data sets for the IT-Model. For the 2T-Model the plot 
are dashed-dotted line (parametrization) with full squares (fitted data sets) corresponding to 
and dashed line (parametrization) with open circles (fitted data sets) to T^. 

V. CONCLUSIONS 



In this work we have suggested a non-equilibrium scenario for the thermal freeze-out in 
heavy-ion collisions. Distinct temperatures are introduced to handle with thermal densities 
of baryons and mesons. Conceptually, these temperatures are associated to the thermal 
state of the vacuum from where particles are produced. In our calculation (Figures [Tjjs]) 
we have shown that within the two-temperatures model we can obtain a rather good fits 
for the ratios of particles population in a wide region of parameters values (Figures (sH). 
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^ (GeV) 

FIG. 5. Plots of baryonic potential against the collision energy. The plot with full line and open 
circles represents the parametrization and fitted data sets for the IT-Model. For the 2r-Model the 
plot with dashed line with open circles is the parametrization and fitted data set. 

A similar situation was also obtained for other analyzed systems. Regarding the fitting of 
particle ratios to the data, we did not privilege any one of the mesonic or baryonic sectors. 
Our results from the two-temperatures model revealed the best fitting when the baryonic 
temperature is higher than the mesonic temperature. Although previous investigations 
about different temperatures for meson and baryon sectors in thermal models have been 
done in regard to Hagedorn's states [28], we remark that in our calculation the attributed 
temperatures are conceptuality differentiated from Hagedorn's temperatures. 

Concerning to the validation of the model, we are aware that further explorations with 
a comprehensive data sets for particles rate are necessary, as well as the analysis of other 
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TABLE VI. Same as Table [Vj this time for the 2r-Model. 



Set 


^SNN (GeV) 


Sm (fm-3) 


Sb (fm-3) 


5*772 1 




^ 1 


CI 


2.32 


0.471 X 10-3 


0.479 X 10-2 


0.0984 


0.167 


1.418 


C2 


17.3 


1.226 


0.524 


2.340 


3.971 


1.465 


C3 


62.4 


1.086 


0.370 


2.937 


3.669 


0.922 


C4 


130 


1.023 


0.379 


2.702 


3.600 


0.929 


C5 


200 


0.974 


0.408 


2.384 


3.568 


0.981 



observables extracted from the experiments. In general, for a given reaction, the set of 
detected particle population includes only a few types of particle ratios and it is a hard task 
to obtain a confident model matching. In addition, it should be noted that it is usual to 
find in the literature different data sets coming from the same experiment, since they are 
constrained to different data filters or selection criteria. Even so, the performed calculation 
(see results in Figures [T]|5| have shown that using the two-temperatures thermal model the 
quality of the fittings for the ratios of particles population were improved in a wide region of 
the parameters values. This opens room enough to accommodate predictions of observable 
not explored here, such as strangeness production, transverse distribution of momentum, 
total particle-antiparticles ratios, and others. 

Finally, we remark that the effect of using the two-temperatures model on other collision 
systems can be enhanced since the conservation laws given by Eqs. (|2])-(|4]) play a central 
role in the calculation. 
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